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Abstract. Astrometric accuracy of complex modern VLBI arrays cannot be calculated analytically. We study 
the astrometric accuracy of phase-referenced VLBI observations for the VLBA, EVN and global VLBI array by 
simulating VLBI data for targets at decimations —25°, 0°, 25°, 50°, 75° and 85°. The systematic error components 
considered in this study are calibrator position, station coordinate, Earth orientation and troposphere parameter 
uncertainties. We provide complete tables of the astrometric accuracies of these arrays for a source separation 
of 1° either along the right ascension axis or along the declination axis. Astrometric accuracy is 50 pas at mid 
declination and is 300 pas at low ( — 25°) and high (85°) declinations for the VLBA and EVN. In extending our 
simulations to source separations of 0.5° and 2°, we establish the formula for the astrometric accuracy of the 
VLBA: AccosS.a = (Aa cos a — 14) x d + 14 (pas) where A^ COS(5 i is the astrometric accuracy for a separation 
d — 1° provided in our tables for various declinations and conditions of the wet troposphere. We argue that this 
formula is also valid for the astrometric accuracy of the EVN and global VLBI array. 

Key words. Astrometry - Techniques: high angular resolution - Techniques: interferometric 



P-i 1. Introduction 



Very Long Baseline Interferometry (VLBI) narrow-angle 
astrometry pioneered by Shapiro et al. (1979) makes use 
of observations of pairs of angularly close sources to cancel 
atmospheric phase fluctuations between the two close lines 
of sight. In this initial approach, the relative coordinates 
between the two strong quasars 3C345 and NRA0512 
and other ancilliary parameters were adjusted by a least- 
squares fit of the differenced phases after connecting the 
VLBI phase s for both sources ov e r a multi-ho ur experi- 
ment. Then, iMarcaide fc Shapiro! (|l983l Il984[) made the 
first phase-referenced map where structure and astrome- 
try were disentangled for the double quasar 1038+528 A 
and B. Both of these experiments demonstrated formal er- 
rors at the level of a few tens of microarcseconds or less in 
the relative angular separation between the two sources. 

Another approach was designed to tackle faint target 
sources by observing a strong reference source (quasar) 
to increase the integ ration time of V LBI from a few min- 
utes to a few hours l|Lestrade et al.lll990|) . This approach 
improves the sensitivity by the factor 
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where Ni, is the number of VLBI baselines, Ti nt is the 
extended integration time permitted by phase-referencing 
(several hours) and T scan is the individual scan length (a 
few minutes). As this factor is very large (e.g. > 50 for 
the 45 baselines of the Very Long Baseline Array), faint 
target sources can be detected and their positions can be 
concomitantly measured with high precision. 

In the approach above, the VLBI phases of the strong 
reference source are connected, interpolated in time and 
differenced with the VLBI phases of the faint source that 
do not need to be connected. The differenced visibilities 
are then inverted to produce the map of the brightness dis- 
tribution of the faint target source and its position is deter- 
mined by reading directly the coordinates of the map peak 
which are relative to the a priori reference source coordi- 
nates. The map is usually highly undersampled but suffices 
for astrometry. This mapping astro metry technique is im- 
plemented in the SPRINT software l|Lestrade et al.l fl990) 
and a similar procedure is also used within the NRAO 
AIPS package to produce phase-referenced VLBI maps 
with absolute source coordinates on the sky. 

While phase-referencing in this way is efficient, it still 
provides no direct positional uncer tainty as does least - 
squares fitting of differenced phases l|Shapiro et alJll979^ . 
In order to circumvent this problem, we have developed 
simulations to evaluate the impact of systematic errors 
in the derived astrometric results. Such simulations have 
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been carried out for of a pair of sources observed with 
the Very Long Baseline Array (VLBA) and the European 
VLBI Network (EVN) at various declinations and angu- 
lar separations. Systematic errors in station coordinates, 
Earth rotation parameters, reference source coordinates 
and tropospheric zenith delays were studied in turn. The 
results of the simulations are summarized below in ta- 
bles that indicate positional uncertainties when consider- 
ing these systematic errors either separately or altogether. 
Such tables can be further interpolated to determine the 
accuracy of any full-track experiment with the VLBA and 
EVN. 

Our study includes atmospheric fluctuations caused by 
the turbulent atmosphere above all stations. These fluc- 
tuations have been considered uniform and equivalent to 
a delay rate noise of 0.1 ps/s for all stations. The impact 
of these fluctuations is limited if the antenna switching 
cycle between the two sources is fast enough. The phase 
structure function measure d at 22 GHz above the VLA by 
ICarilli Holdawavl Jl999) provides prescriptions on this 
switching tim e. At high f r equen cy, it can be as short as 
10s, as e.g. in iReid et alJ l|2003|) who carried out precise 
43 GHz VLBA astrometric observations of Sgr A* at a 
declination of —28°. Switching time in more clement con- 
ditions is typically a few minutes at 8.4 GHz for northern 
sources. 

A few applications of mapping astrometry are the 
search for extra-solar planets around radio-emitting stars 
(jLestrade et al.lll994l) . the determin ation of the Gravit y 
Probe B guide star proper motion ijLebach et alJll999l) . 
the determination of absolute motions of VLBI compo- 
nents i n extragalactic sour ces, e.g. in compac t symetric 
objec ts ijCharlot et al . 2005) or core-jet sources l|Ros et alJ 
Il999|) . probing the jet collimation region in extragalactic 
nuclei l)Lv et all l2004h . pulsar parallax and proper mo- 
tion measurements l(Brisken et al.l Eo02'l and the determi- 
nation of parallaxes and proper motions of maser sources 
in the whole Galaxy as planned with the V ERA project 
l|Kawaguchi et alJl200d: iHonma et aLll2000j) . 



systematic errors in the model of the phase limit narrow- 
angle astro metry precision to roug hly ten times this level 
in practice ijFomalont et all Il999) . An analytical study 
of systematic errors in phase-refere nced VLBI astrome- 
try over a single baseline is given in IShapiro et all l|l979t) 
and it shows that all systematic errors are scaled by the 
source separation. Another error analysi s in such differen - 
tial VLBI measurements can be found in lMorabitol l)l984|) . 
However, for modern VLBI arrays with 10 or more anten- 
nae, the complex geometry makes the analytical approach 
intractable. For this reason, we have estimated such sys- 
tematic errors by simulating VLBI visibilities and invert- 
ing them for a range of model parameters (station coor- 
dinates, reference source coordinates, Earth Orientation 
parameters, and tropospheric dry and wet zenith delays) 
corresponding to the expected errors in these parameters. 

The visibilities were simulated for a pair of sources at 
declinations -25°, 0°, 25°, 50°, 75°, 85° and with angu- 
lar separations 0.5°, 1° and 2° for the VLBA, EVN and 
global VLBI array (VLBA+EVN). For each of these cases, 
we simulated visibilities every 2.5 min from source rise to 
set (full track) with a lower limit on elevation of 7°. The 
adopted flux for each source (calibrator and target) was 
1 Jy to make the phase thermal noise negligeable in our 
simulations. For applications to faint target sources, one 
should combine the corresponding thermal astrometric un- 
certainty (Eq. 1) with the systematic errors derived below. 
The simulated visibilities were then inverted using uniform 
weighting to produce a phase-referenced map of the tar- 
get source and estimate its position. This operation was 
repeated 100 times in a Monte Carlo analysis after vary- 
ing slightly the parameters of the model based on errors 
drawn from a Gaussian distribution with zero-mean and 
plausible standard deviation. We report the rms of the 
differences found between the known a priori position of 
the target source and the resulting estimated positions as 
a measure of the corresponding systematic errors for each 
of the above cases. We have adopted the usual astrometric 
frequency of 8.4 GHz for this analysis. 



2. Method 

As indicated in e.g. lThomnson et al.l 1 1986). the theoreti- 
cal precision of astrometry with the interferometer phase 
is 

1 1_ A_ 

2tt SNR B' 



&a,5 = 7T~ 



(1) 



where SNR is the signal-to-noise ratio of the observa- 
tion, A is the wavelength and B is the baseline length 
projected on the sky. For observations with the VLBA 
(B ~ 8000 km), A = 3.6 cm, and a modest SNR of 10, this 
theoretical precision is breathtakingly ~ 15 /xas. Although 
a single observation of the target yields an ambiguous po- 
sition, multiple observations over several hours easily re- 
move ambiguities eve n with a sparse u-v plane coverage 
l|Lestrade et al.lll990j) . 

While the theoretical precision above might be re- 
garded as the potential accuracy attainable for the VLBI, 



3. Phase model used in simulation 

T he phase delay and group delay in VLBI are described 
in lSovers et al.l l|l998l) . The phase <fi = vt at frequency v 
is related to the interferometer delay 



Ttr 



TR+T s tr 



Talk- 



Specifically, the geometric delay is: 
r g = [P][N}[EOP] 



b.k 



with the precession matrix [P], the nutation matrix [N], 
the Earth Orientation Parameters matrix [EOP], the 
baseline coordinates b in the terrestrial frame, the source 
direction coordinates k computed with source right ascen- 
sion and declination in the celestial frame. The "retarded 
baseline correction" to account for Earth rotation dur- 
ing elapsed time r g must also be modelled l(Sovers et alJ 
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Table 1. Adopted rms errors for the source coordi- 
nates, VLBA station coordinates and Earth Orientation 
Parameters in our Monte Carlo simulations. 



Parameters 


Errors 


Source coordinates 




ao cos So 


0.25/1 mas 


do 


0.25/1 mas 


Station coordinates 




X 


1-2 mm 


Y 


1-3 mm 


Z 


1-2 mm 


Earth Orientation Parameters 




Ap 


0.2 mas 


Y P 


0.2 mas 


UT1 - UTC 


0.02 ms 


ip sin e 


0.3 mas 


e 


0.3 mas 



Il998l) . The differential tropospheric delay T trop between 
the two stations is computed with a static tropospheric 
model and the simple mapping function 1 / sin E (where E 
is the source elevation at station) to transform the zenith 
delay into the line-of-sight delay at each station. The dif- 
ferential ionospheric phase delay T iono — — fcTEC/^ 2 is 
related to the total electronic content TEC in the direc- 
tion of the source at each station. The General Relativity 
delay tr takes into account light propagation travel time 
in the gravitational potential of the Sun. The source struc- 
ture cont ribution Tw r ,,, r . can be computed according to the 
model by [Chariot] l)l99d(l but was not included in our sim- 
ulations which are for point sources. The clock delay r c ik 
cancels in differenced VLBI phases. 

The model above is that implemented in the SPRINT 
software used for our simulations. It is thought to be com- 
plete for narrow-angle astrometry and additional refine- 
ments, such as ocean loading, atmospheric loading, etc. . . , 
would not make difference into our results. We have not 
studied the ionosphere contribution to systematic errors. 
The unpredictible nature of the ionosphere makes this task 
difficult. Calibration of the ionosphere by dual- frequency 
observations, or ove r a wide bandwidth at low frequency 
^Brisken et al.ll2002|) . or simply by observing at high fre- 
quency (> 10 GHz) where the effect is small, offers solu- 
tions to this problem. 

4. Results 

4.1. VLBA 



The parameter rms errors adopted as plausible for the 
VLBA phase model are listed in Tables Q and The ref- 
erence source coordinate uncertainties (Aaocos^o, A5q) 
of 1 ma s are typical of tho se in the VLBA Calibrator 
burvey ()Beaslev et alJ 12002^1 . from which most of the 
reference sources originate. However, ICRF extragalactic 
sources have be tter position accuracies down to 0.25 mas 
I Ma et al.ll998l) . We have thus carried out the calculations 
for both of these cases (1 mas and 0.25 mas) and both a 



Table 2. Dry and wet tropospheric zenith path delays 
(i~dtrp and T wtr p) at the VLBA stations along with the 
adopted rms errors Ar^trp and AT wtrp in our Monte Carlo 
simulations. 



Stations 



Dry trop. 



Wet trop. 
Mean Max 





7~dtrp 




'Twtrp 


^Twtrp 


Twtrp 


^Tyjtrp 




(cm) 


(cm) 


(cm) 


(cm) 


(cm) 


(cm) 


Brewster 


225 


0.5 


8 


2.7 


13 


4.3 


Fort Davis 


192 


0.5 


8 


2.7 


15 


5.0 


Hancock 


223 


0.5 


9 


3.0 


19 


6.3 


Kit Peak 


185 


0.5 


6 


2.0 


15 


5.0 


Los Alamos 


185 


0.5 


6 


2.0 


13 


4.3 


Mauna Kea 


149 


0.5 


1 


2.0 


4 


2.0 


North Liberty 


225 


0.5 


10 


3.3 


19 


6.3 


Owens Valley 


199 


0.5 


5 


2.0 


20 


6.7 


Pietown 


176 


0.5 


4 


2.0 


12 


4.0 


Saint Croix 


213 


0.5 


22 


7.3 


30 


10.0 



and <5o have been perturbed by these uncertainties in our 
simulations. The uncertainties for the station coord inates 
are from the ITRF2000 frame ^Boucher et al.ll2004h while 
those for the Earth Orientation Parameters are from the 
IERS web site 1 . The adopted dry tropospheric rms error 
A.Tdtrp of 0.5 cm corresponds to 2.5 millibars in atmo- 
spheric pressure uncertainty at sea level. Although barom- 
eter reading is usually better, the absolute calibration of 
station barometers is at this level. Uncertainties in the 
wet tropospheric zenith delay T w t rp derive d from temper- 
ature and humidity are known to be large ( SaastamoinerJ 
Il973j) . Experience makes us believe that a 30% error is 
likely on T wtrp and thus we took 1/3 of T w t rp as the plau- 



sible rms error Ar w t rp with a minimum value of 2 cm. We 
carried out simulations for both mean and maximum val- 
ues of wet zenith path delays based on estimates of T wtrp 
recently derived from multipl e VLBA geodetic and astro- 
metric sessions llSover.jbon^ . The maximum wet zenith 
delays and corresponding errors were used to investigate 
the impact of extreme weather conditions on observations. 
These values are listed for each VLBA station in Tabled 
We simulated the visibilities of a full u-v track ex- 
periment with the VLBA for six declinations between 
—25° and 85° with a 1° relative source separation (ei- 
ther oriented in right ascension or in declination) . Uniform 
weighting was applied to the visibilities, resulting in a syn- 
thesized beam mainly shaped by the longest baselines. 
As a test, we have also removed the 9 baselines smaller 
than 1500 km in length out of the 45 baseline array and 
noted a decrease in systematic errors of ~ 15% in a few 
test cases. Conservatively, we have retained these "short" 
baselines in our final simulations. This is motivated by the 
fact that all possible baselines must be kept for sensitiv- 
ity when observing weak sources. The antenna switching 
cycle between target and reference sources was set to 2.5 



1 http: / /hpiers.obspm.fr/iers /eop /eopc04/EOPC04. GUIDE 
(Table 2). 
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Table 3. VLBA rms astrometric errors (in /zas) for a relative source separation (a — ao)cos<5o = 1°. Individual 
astrometric error contributions from calibrator position, Earth orientation parameter, station coordinate, and dry and 
wet troposphere uncertainties are given separately, while the last two lines indicate the total astrometric errors when 
all model parameters are perturbed together. 



Declination of source 





-25° 




0° 




25° 




50° 




75° 




85° 




Error component 


Aq cos S 


AS 


Aq cos S 


AS 


Aa cos S 


AS 


Aa cos S 


AS 


Aa cos S 


A<5 


Aa cos S 


AS 


Calibrator position (1 mas error) 


8 


7 


1 


9 


8 


16 


20 


26 


59 


68 


196 


193 


Calibrator position (0.25 mas error) 


2 


7 


1 


3 


2 


5 


3 


5 


12 


11 


49 


50 


Earth orientation 


1 


8 


1 


5 


1 


6 


1 


5 


1 


4 


1 


4 


Antenna position 


2 


8 


2 


4 


2 


4 


2 


3 


2 


3 


2 


3 


Dry troposphere 


15 


45 


9 


16 


7 


9 


10 


11 


18 


23 


14 


16 


Wet troposphere (mean) 


53 


182 


34 


57 


33 


28 


31 


45 


54 


72 


79 


88 


Wet troposphere (max) 


87 


219 


46 


66 


42 


38 


49 


56 


65 


78 


81 


91 


Total (mean wtrp) 


60 


175 


36 


50 


33 


32 


37 


53 


87 


103 


227 


258 


Total (max wtrp) 


85 


217 


43 


74 


42 


44 


46 


66 


100 


117 


226 


240 



Table 4. VLBA rms astrometric errors (in fxas) for a relative source separation S — So = 1°. Individual astrometric 
error contributions from calibrator position, Earth orientation parameter, station coordinate, and dry and wet tropo- 
sphere uncertainties are given separately, while the last two lines indicate the total astrometric errors when all model 
parameters are perturbed together. 



Declination of source 





-25° 




0° 




25° 




50° 




75° 




85° 




Error component 


Aa cos S 


AS 


Aa cos S 


AS 


Aa cos S 


AS 


Aa cos S 


AS 


Aa cos S 


AS 


Aa cos S 


AS 


Calibrator position (1 mas error) 


7 


8 


1 


7 


8 


11 


21 


2 


59 


2 


199 


1 


Calibrator position (0.25 mas error) 


2 


7 


1 


3 


2 


4 


5 


2 


20 


1 


43 


1 


Earth orientation 


5 


7 


5 


3 


5 


4 


4 


3 


3 


1 


3 


1 


Antenna position 


2 


9 


2 


6 


2 


5 


3 


3 


2 


2 


2 


2 


Dry troposphere 


17 


54 


6 


19 


2 


12 


5 


9 


11 


9 


12 


13 


Wet troposphere (mean) 


80 


272 


32 


98 


11 


41 


12 


32 


43 


41 


60 


62 


Wet troposphere (max) 


112 


358 


43 


114 


19 


61 


17 


46 


59 


55 


74 


71 


Total (mean wtrp) 


84 


284 


30 


99 


16 


42 


25 


33 


81 


36 


189 


67 


Total (max wtrp) 


121 


481 


44 


134 


20 


56 


26 


46 


92 


65 


212 


74 



minutes. The results, however, do not depend critically 
on this value. It was chosen so that the automatic phase 
connection routine for the reference source does not dis- 
card too much data in the presence of a delay rate error 
of 0.1 ps/s (adopted uniformly for all the stations in the 
simulation). As mentioned previously, we analysed these 
data simulated with SPRINT using the a priori parame- 
ter values perturbed by some errors. We carried out this 
analysis 100 times for each systematic error component 
with perturbation errors drawn from Gaussian distribu- 
tions with zero mean and standard deviations according 
to the rms errors in Tables QJand [21 The resulting position 
of the target was estimated by reading the peak position 
in each of the 100 phase-referenced maps. The pixel size 
in the maps was 0.05 mas. This size is small compared to 
the synthesized beam (~ 1 mas at 8.4 GHz on 8000 km 
baseline) and, hence, the uncertainty in the peak posi- 
tion due to the pixel size is negligeable. This position was 
determined by fitting a parabola over the full half beam 
width. This procedure was used in the Hipparcos/VLBI 



work of iLestrade et all l|l999|) and was found to be ap- 
propriate. As expected, each position was slightly offset 
from the map phase center, reflecting the corresponding 
systematic errors. After substracting the initial pertur- 
bation in the calibrator position, we calculated the rms 
of these 100 relative coordinate offsets AacosS and AS 
for the adopted 1° source separation in right ascension or 
declination. Note that the mean of these 100 coordinate 
offsets was close to zero in all cases. In Tables and 0] 
we report the rms astrometric errors for each individual 
error component along with the total astrometric errors 
when all model parameters are perturbed together in the 
simulation. The total errors were derived by considering a 
1 mas error in the calibration position. 

The wet troposphere systematic error clearly domi- 
nates over all the other error components for S < 50° but 
the calibrator error dominates at higher declinations if its 
position is not kno wn to better than 1 m as. This behavior 
was first noted bv IShaniro et al.l l|l979F ) who derived an- 
alytical formulae providing the astrometric errors caused 
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Table 5. Reduced chi-square xl an d probability p of 
Gaussian distribution for the astrometric errors Aa cos S 
and AS using the Pearson test. 







(a — OQ 


cos So = 


1° 




Declination 




Aa cos 5 






AS 






xl 


P 




xl 


P 




-25° 


1.51 


1.43x10^ 




2.10 


9.59x10" 


-6 


0° 


1.67 


1.29xl0" 2 




1.73 


1.95x10" 


-3 


25° 


2.33 


8.12xl0" 8 




1.97 


1.12x10 


-3 


50° 


0.97 


0.511 




1.18 


0.309 




75° 


0.84 


0.783 




1.22 


0.273 




85° 


1.06 


0.659 




1.02 


0.424 








(5- 


So) = 1° 






Declination 




Aa cos 8 






AS 






xl P 


xl 


P 




-25° 


1.19 


0.286 




0.76 


0.885 




0° 


1.62 


2.24xl0" 2 




1.35 


6.01x10 


-2 


25° 


3.22 


4.62xl0" 4 




1.91 


5.57x10" 


-4 


50° 


2.23 


5.19xl0~ 6 




2.18 


1.28x10 


-4 


75° 


1.38 


0.116 




2.05 


1.99x10" 


-3 


85° 


1.22 


0.260 




2.98 


5.62x10" 
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by the calibrator coordinate uncertainties in the case of a 
single VLBI baseline. A detailed analysis comparing our 
simulated errors with those obtained from these formulae 
is given in Appendix A. Other systematic errors, in partic- 
ular the Earth orientation parameter and the station co- 
ordinate errors, are small. In Tables [3] and we note that 
astrometric errors originating from mean and maximum 
wet troposphere uncertainties are not drastically different 
(a ratio of 1.5 at most). 

Finally, we have plotted in Fig. Q] the distribution of 
all coordinate offsets Aa cos 5 and AS for the 50° declina- 
tion target when all perturbation errors are present. For 
this specific case we have carried out 1000 simulations to 
refine the binning of the distribution. We have also per- 
formed the Pearson test on all distributions and provide 
the reduced chi-square xl and probability p that such dis- 
tributions are Gaussian in Tabled The results of this test 
show that most of the distributions are not Gaussian with 
p generally smaller than 0.4. 

4.2. EVN 

We have carried out a similar study for the EVN by simu- 
lating full track observations for the 10 stations of the ar- 
ray at 8.4 GHz. The adopted errors for the reference source 
coordinates and Earth orientation parameters were iden- 
tical to those used in the VLBA simulations. Station co- 
ordinate errors were similar to the VLBA ones (1-6 mm), 
with the exception of those for Weste rbork which are at 
the level of 50 mm [|Charlot et all2002|) . The same scheme 
as that adopted for the VLBA was used to define zenith 
dry and wet tropospheric delay errors at each EVN station 
and the corresponding values are given in Table El 

Since the EVN comprises antennas with different sen- 
sitivities, each baseline has been weighted by the recip- 



Table 6. Dry and wet tropospheric zenith path delays 
(Tdtrp and T wtrp ) at the EVN stations along with the 
adopted rms errors Ar^trp and AT wtrp in our Monte Carlo 
simulations. 



Stations Dry trop. Wet trop. 

Mean Max 





7~dtrp 


ATdtrp 


7~wtrp 


^Twtrp 


'Twtrp 


^Twtrp 




(cm) 


(cm) 


(cm) 


(cm) 


(cm) 


(cm) 


Effelsberg 


220 


0.5 


8 


2.7 


20 


6.7 


Hartebeesthoek 


199 


0.5 


10 


3.3 


17 


5.7 


Medicina 


231 


0.5 


11 


3.7 


18 


6.0 


Noto 


229 


0.5 


12 


4.0 


20 


6.7 


Onsala 


230 


0.5 


8 


2.7 


14 


4.7 


Sheshan 


231 


0.5 


22 


7.3 


36 


12.0 


Urumqi 


210 


0.5 


10 


3.3 


10 


3.3 


Westerbork 


220 


0.5 


8 


2.7 


20 


6.7 


Wettzell 


215 


0.5 


7 


2.3 


13 


4.3 


Yebes 


208 


0.5 


5 


2.0 


5 


2.0 



rocal of their noise power equivalent VSEFDi x SEFD 2 
with System Equivalent Flux Densities (SEFD^) for each 
station according to Table 2 of the EVN Status Table 2 
(as available in May 2003). The Effelsberg-Westerbork 
baseline is the most sensitive baseline of the array but 
also the shortest one and so unfavorable for high-accuracy 
astrometry. For this reason, we decided to perform the 
simulations without this baseline, hence using an array 
of 44 baselines only. We have applied uniform weighting 
to the visibilities similarly to the VLBA. We have tested 
that in removing the 12 baselines shorter than 1500 km 
in this 44 baseline array, systematic errors decrease by 
~ 20% but, conservatively, we have kept them in our sim- 
ulations. In order to reduce the number of simulations, 
calculations were carried out for only mean values of the 
wet zenith tropospheric delays since the results when us- 
ing mean or maximum values were not found to be drasti- 
caly different. We also did not calculate individual contri- 
butions from calibrator position, dry tropospheric zenith 
delay and Earth orientation parameter errors since these 
were found to be very small for the VLBA (sec Tables |31 
and^l. One should keep in mind, however, that calibra- 
tor error dominates at high declination. The results of the 
EVN simulations are reported in Tables [7| and [S] for a 1° 
source separation in right ascension or declination. 

At declination —25°, many SPRINT maps were found 
to be ambiguous, i.e. the main lobe of the point spread 
function of the EVN could not be identified because sec- 
ondary lobes were too high. This is essentially caused by 
the relatively high latitude of the array and hence to the 
difficulty of observing such low declination sources due to 
very limited visibility periods. For this reason, we do not 
provide EVN results for this declination. For other de- 
clinations, EVN astrometric errors (Tables and [SJ are 
similar to those found for the VLBA (Tables |3] and 
and the Westerbork position error is not a limiting factor. 

2 http://www.mpifr-bonn.mpg.de/EVN/EVNstatus.txt 
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-100 100 

Acx*cos(d) (^as) A6 (/ids) 

Fig. 1. Distribution of total astrometric errors for a 1° relative source separation along declination at declination 50°. 
All perturbating errors (calibrator position, Earth orientation parameters, station coordinates, dry and wet troposphere 
parameters) are considered together in this simulation. 

Table 7. EVN rms astrometric errors (in /^as) for a relative source separation (a — ao)cos5o = 1°- Individual 
astrometric error contributions from station coordinate and wet troposphere uncertainties are given separately along 
with the total astrometric errors when all model parameters are perturbed together. 











Declination of source 












0° 




25° 




50° 


75° 




85° 




Error component 


Aa cos 8 


AS 


Aa cos 8 


AS 


Aa cos 8 AS 


Aa cos S 


AS 


Aa cos 8 


AS 


Antenna position 


5 


4 


5 


4 


6 5 


5 


5 


5 


5 


Wet troposphere (mean) 


55 


11 


37 


14 


52 33 


73 


40 


65 


31 


Total (mean wtrp) 


57 


12 


44 


15 


57 45 


91 


81 


206 


185 



Table 8. EVN rms astrometric errors (in fias) for a relative source separation 6 — 5q = 1°. Individual astrometric 
error contributions from station coordinate and wet troposphere uncertainties are given separately along with the total 
astrometric errors when all model parameters are perturbed together. 











Declination of 


source 












0° 




25° 




50° 




75° 




85° 




Error component 


Aa cos 8 


AS 


Aa cos S 


A<5 


Aa cos 5 


AS 


Aa cos S 


AS 


Aa cos S 


AS 


Antenna position 


7 


6 


7 


5 


4 


5 


5 


5 


4 


5 


Wet troposphere (mean) 


51 


31 


29 


54 


18 


81 


31 


58 


33 


68 


Total (mean wtrp) 


62 


29 


33 


57 


27 


78 


79 


61 


201 


61 




Declination accuracies are somewhat better for the EVN 
than for the VLBA at low declination (0° and 25°), a 
consequence of the participation of Hartebeeshoek (South 
Africa) in such observations. 

4.3. Global VLBI array 

We have carried out a similar study for the global VLBI 
array which is the combination of the VLBA and EVN. 
It includes 20 stations, with 190 possible baselines. As 



discussed above, the Effelsberg-Westerbork baseline was 
ignored and the calculations were thus carried out for 
189 baselines only. The adopted systematic error values 
for the simulations with this array were the same as those 
adopted for the individual VLBA and EVN (Tables □ 
and|HJ) and calculations were performed for full track ob- 
servations as previously. The results of these simulations 
( Tables l9l and ITU)> indicate that the astrometric errors for 
the global VLBI array are consistent with those found for 
the VLBA and the EVN. As expected, these errors are 
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Table 9. Global VLBI array rms astrometric errors (in /xas) for a relative source separation (a — ao) cos<5o = 1°. The 
individual astrometric error contribution from wet troposphere uncertainties is given separately along with the total 
astrometric errors when all model parameters are perturbed together. 













Declination of 


source 














-25° 




0° 




25° 




50° 




75° 




85° 




Error component 


Aa cos S 


A<5 


Aa cos S 


A<5 


Aa cos 5 


AS 


Aa cos S 


A<5 


Aa cos S 


A<5 


Aa cos S 


A<5 


Wet troposphere (mean) 


71 


76 


32 


42 


26 


34 


23 


13 


22 


6 


27 


9 


Total (mean wtrp) 


82 


67 


34 


46 


24 


44 


34 


33 


64 


76 


196 


203 



Table 10. Global VLBI array rms astrometric errors (in fxas) for a relative source separation 6 — Sq = 1°. The 
individual astrometric error contribution from wet troposphere uncertainties is given separately along with the total 
astrometric errors when all model parameters are perturbed together. 













Declination of source 














-25° 




0° 




25° 




50° 




75° 




85° 




Error component 


Aa cos S 


AS 


Aa cos 5 


AS 


Aa cos S 


AS 


Aa cos S 


AS 


Aa cos 8 


AS 


Aa cos S 


AS 


Wet troposphere (mean) 


60 


305 


26 


71 


10 


43 


9 


44 


5 


21 


5 


26 


Total (mean wtrp) 


61 


279 


24 


78 


15 


45 


22 


46 


61 


24 


183 


27 



generally slightly better than the ones derived for each 
individual array. 

5. Discussion 

5.1. General results 

Our simulations show that the astrometric accuracy 
of the VLBI phase-referencing technique (defined as 
y/ (Aa cos S) 2 + (AS) 2 ) is ~ 50 /ias for mid declinations 
and is < 300 /j,as at low and high declinations for point 
sources with a relative separation of 1°. The major system- 
atic error components are the wet tropospheric delay and 
the calibrator astrometric position, the latter only at high 
declination. Station coordinate, Earth orientation param- 
eter and dry tropospheric zenith delay errors contribute 
generally to less than 20 /xas in the error budget. 

5.2. Simulation of the VLBA without Saint Croix 

We speculated that if the VLBA station at Saint Croix in 
the Virgin Islands that suffers from dampness were with- 
drawn from the array, it should improve the astrometric 
accuracy of the VLBA. We thus repeated our VLBA sim- 
ulations without that station. The results of this test are 
given in TablcsEl an dEl In contrast to our intuition, the 
astrometric accuracy is actually degraded when the target- 
calibrator direction is oriented along declination. In fact, 
the addition of Saint Croix strengthens the geometry of 
the array and improves the astrometric accuracy despite 
severe weather conditions. In order to further explore this 
question, we ran simulations without Pic Town in the mid- 
dle of the array and without Mauna Kea at the far West 
of the array. Withdrawing Pic Town does not change the 
astrometric accuracy but the absence of Mauna Kea de- 
grades the accuracy in a similar way to Saint Croix. 



5.3. Linearity of the astrometric accuracy with source 
separation 

An important question is whether the astrometric accu- 
racy scales linearly as a function of the source separation. 
To study this matter, we repeated all the previous simu- 
lations but with source separations of 0.5° and 2°. Then, 
we performed a linear fit to the astrometric errors for the 
three values of the calibrator-target separation (0.5°, 1° 
and 2°), considering separately each systematic error com- 
ponent of the tables above. Figure shows an example of 
such results for the VLBA in the case of a target at +25° 
declination. Overall, our plots show that the astrometric 
accuracy generally scales fairly linearly as a function of 
the source separation. 

To obtain a quantitative measure of the likchood of the 
linearity, we determined the regression coefficients for each 
of the 107 linear fits. Such coefficients should be close to 
1 for a linear behavior while they should decrease as the 
behavior becomes less linear. This analysis reveals that 
80% of the coefficients are larger than 0.95, indicating that 
the astrometric errors behave linearly. Among all errors, 
calibrator position systematics are those that were found 
to behave the least linearly. An empirical formula for the 
astrometric accuracy A acos s,8 has been further estimated 
by averaging the parameters of all the fits : 

A a cos 5 , 5 = ( Al° cos Sfi - 14) x d + 14 (fxas) (2) 

where A* cosi5 s is the astrometric error for 1° source sep- 
aration as provided by our tables J31and0]for the VLBA, 
EandElfo r the EVN, ED and ITO! for the gl obal VLBI array) 
and d — yj ((a — ao) cos<5 ) 2 + (S — So) 2 is the source sep- 
aration in degrees. In Section l4~2l we noted that the astro- 
metric accuracies of the EVN and the VLBA are similar, 
hence this formula should apply to the EVN, too. 



8 



Pradel, Chariot, Lestrade: Astrometric accuracy of phase-referenced observations with the VLBA and EVN 



Table 11. VLBA without Saint Croix rms astrometric errors (in /xas) for a relative source separation (a — ao) cos Jo = 
1°. The individual astrometric error contribution from wet troposphere uncertainties is given separately along with 
the total astrometric errors when all model parameters are perturbed together. 













Declination 


of 


source 














-25° 




0° 




25° 






50° 




75° 




85° 




Error component 


Aa cos S 


AS 


Aa cos S 


AS 


Aa cos 5 


AS 




Aa cos S 


A<5 


Aa cos S 


A<5 


Aa cos 5 


A<5 


Wet troposphere (mean) 


62 


171 


28 


47 


27 


57 




34 


65 


46 


63 


37 


67 


Total (mean wtrp) 


63 


193 


27 


41 


31 


62 




42 


68 


83 


82 


211 


207 



Table 12. VLBA without Saint Croix rms astrometric errors (in /xas) for a relative source separation (5 — <5o = 1°. The 
individual astrometric error contribution from wet troposphere uncertainties is given separately along with the total 
astrometric errors when all model parameters are perturbed together. 













Declination of 


source 














-25° 




0° 




25° 




50° 




75° 




85° 




Error component 


Aa cos S 


AS 


Aa cos 5 


AS 


Aa cos S 


AS 


Aa cos S 


AS 


Aa cos S 


AS 


Aa cos 8 


AS 


Wet troposphere (mean) 


183 


563 


62 


189 


19 


55 


17 


41 


39 


27 


42 


44 


Total (mean wtrp) 


190 


534 


70 


191 


24 


71 


26 


42 


74 


28 


216 


40 



45 
40 

35 - 

8 30 - 



25 



20 



Calibrator (1 mas) — I— 

Calibrator (0.25mas) --X-- 
EOP 

Stations -S- 

Dry troposphere -#- 

Wet troposphere -~0~ 
All 




0.5 1 

Relative source separation ( 



60 



10 



Calibrator (1 mas; 

Calibrator (0.25mas) ~-X 
EOP 

Stations □ 

Dry troposphere --■ 

Wet troposphere -~Q 

All • 




0.5 1 1.5 

Relative source separation f) 



Fig. 2. Astrometric accuracy as a function of the relative source separation for a target-calibrator pair observed 
with the VLBA at a declination of +25°. Each error component is represented with a different symbol and the total 
contributing error is also plotted. 



As a verification of this empirical formula, we com- 
puted the astrometric accuracy for eight target-calibrator 
pairs observed with the global VLBI array as part of a 
project to moni tor absolute lobe mo tions in compact sym- 
metric objects dCharlot et alJ 120051) . For the source pair 
J2212+0152/J2217+0220 with a separation of 1.37° along 
the right ascension direction, we obtained simulated ac- 
curacies Aa cos 5q = 42 fias and Ad = 63 /xas, versus 
Aa cos 5q = 44 fias and AS = 63 fias when derived from 
Eq. [3 and Table El In the worst case (target-calibrator 
J0754+5324/J0753+5352 with a separation of 0.50° along 



declination), simulated accuracies were Aa cos Sq — 18 /ias 
and AS = 12 ^as while Eq. and Table E3 give 
AacosJo = 20 fia.s and A5 = 24 /ias. Thus, overall we 
found a discrepancy of a factor of 2 at most between our 
simple formula (Eq. |2J) and real simulation of the case 
considered. 

6. Conclusion 

We have performed extensive simulations of VLBI data 
with the VLBA, EVN and global VLBI array to study 
the dependence of the astrometric accuracy on systematic 
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errors in the phase model of phase-referenced VLBI ob- 
servations. Systematic errors considered in this study are 
calibrator position uncertainties, station coordinate uncer- 
tainties, Earth Orientation Parameters uncertainties and 
dry and wet troposphere errors. We have adopted state of 
the art VLBI values for these errors. 

Our simulations show that the astrometric accuracy of 
a full track phase-referenced VLBI experiment is 50 /xas 
at mid declination and is ~ 300 /xas at low (—25°) and 
high (+85°) declinations for point sources angularly sepa- 
rated by 1°. Not surprinsingly, the major systematic error 
originates from wet tropospheric zenith delay uncertain- 
ties except at high declination where calibrator position 
uncertainties dominate. We show that the astrometric ac- 
curacy A acos g s depends linearly on the source separa- 
tion and we established the simple formula A a cos s.s — 
( A aco S <5,5 - 14) x d+14 (fias) where A^ cosS>s is the astro- 
metric error provided by our tables for the various arrays 
and configurations and d — y/ ((a — ao) cos S) 2 + (S — So) 2 
is the source separation in degrees. Our study has been 
carried out for point sources but variable source structure 
is likely to degrade the accuracy derived from this formula. 



Appendix A: Analytical Behavior 

The analytic al formulae in the Appendix A of 

IShapiro et alJ l)l979|) provide the astrometric errors caused 
by the inaccuracy of the calibrator coordinates in the case 
of a single VLBI baseline. Adopting our notation, these 
formulae become : 



Aa ~ 



and 



((a — ao) tan5)A5o 

-((5-5 )tan5+l/2x (a - a ) 2 )Aa , 



AS ~ (-(5-5 )cot5 + l/2x {a - a Q ) 2 )A8 + 
((a — ao) cot 5) Aao. 

where Aa and A5 are the errors in right ascension and 
declination introduced by errors Aao and A5o in the co- 
ordinates of the reference source. The expression above for 
A6 restores correctly the last term of the equation which 
was misprinted in the original paper. These simple for- 
mulae are, however, valid only for the special geometry 
adopted by the authors where the "baseline declination" 

o 

is . 

Adopting the same parameters as in our simulations 
(Aao = l/cos5o mas, ASq = 1 mas, a — ao = 0° or 
(1/ cos<5o)°, <5 — (5o = 1° or 0°), we obtain the astromet- 
ric errors plotted as a function of declination in Fig IA.1I 
(dotted lines). The results of our simulations for declina- 
tions of -25°, 0°, 25°, 50°, 75° and 85° in the case of the 
VLBA (first lines of Tables 3 and 4) are also superimposed 
on these plots. 

The right ascension errors obtained from the simula- 
tions match perfectly those derived analytically, while the 
declination errors show a strong discrepancy near declina- 
tion 0° (although they agree at high declinations). This 
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Fig. A.l. Astrometric errors Aacos5 and A5 (respec- 
tively left and right) as a function of declination. The 
two upper plots are for the case a — a = (l/cos5)° and 
5 — So = 0° while the two lower plots are for the case 
a — a = 0° and S — do = 1°. T he continuous dotted lines 
show the errors derived from the lShapiro et all l|l979l) for- 
mulae. The stars show the errors from our simulations at 
six declinations from —25° to 85°. 



discrepancy originates from a singularity in the AS for- 
mula at S = 0° (term in cot 5) , inherent to the approxi- 
mation used to establish the formula (baseline declination 
of 0°). For a more complex and realistic network, such 
a singularity does not exist, as also demonstrated by the 
results of our simulations. 
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